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An  I ones  in  many  countries  li.ni-  to  lace  problems  ol  aeroelasiieity  .nul  tint  lot  with 
atrcialt  carry  mg  more  am!  more  stores  Flu*  two  papers  ol  this  Report  presented  to  the  Snh 
Committee  on  Aeroelasiieity  ol  the  Ntrnetutvs  amt  Materials  Panel  during  the  doth  Meeting 
ol  the  Panel  ileal  with  two  different  aspects  ol  the  prohlem 

I he  paper  ol  A l ot/e  is  concerned  with  the  many  configurations,  symmetric  or  asymmetric, 
that  may  occur  and  with  then  .onsei|uences  on  the  natural  modes  ol  the  structure  It  clarifies 
the  difficulties  one  has  to  face  to  clear  the  (light  domain  lor  all  (light  configurations  ami 
proposes  useful  approaches 

1 lie  paper  ol  C Hwang.  B.A  Wmtlier.  1.1  Noll  and  M li  l-'armer  deals  with  the  difficult 
problem  ol  tint  tor  suppression,  it  exposes  the  approach,  the  design,  and  the  wind  tunnel  tests 
ol  the  model  ol  a fighter,  carry  ing  stores  and  equipped  with  a tint  ter  suppression  device 

Ihe  two  papers  are  of  great  inicrvsf  to  aeroelasttctans  and  may  goe  useful  help  to  the 
designer. 
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ASYMMETRIC  STORE  FLUTTER 


by 


A.  Lotze 

MESSERSCHMITT-BOLKOW-BLOHM  GmbH. 
Unternehmenaberelch  Flugzeuge 
Postfach  801160  - 8 MUnchen  80 
W. -Germany 


INTRODUCTION 

Numerous  missions  of  military  aircrafts  are  dealing  with  asymmetrical  store  configu- 
rations. External  stores  may  be  flown  asymmetrically  throughout  the  whole  mission  like 
electronic  pods,  camera  pods  and  pylon  mounted  fuel  tanks  In  combination  with  external 
weapon  carriage.  Asymmetrical  store  configurations  also  can  occur  temporarily  by  asym- 
metrical sequences  of  weapon  release. 

1 Asymmetrical  effects  generated  by  tolerances  In  mass  and  stiffness  distributions 

are  generally  expected  to  be  less  lmi>ortant  for  flutter,  but  may  create  considerable 
problems  during  ground  resonance  and  flight  flutter  tests.  Caused  by  closely  spaced 
f requencles,  beat ing  of  resonance  mode  frequencies  can  arise  which  makes  it  difficult 
to  Identify  the  mode  and  to  evaluate  the  exact  damping. 


Performing  flutter  work  it  is  quite  common  to  base  the  investigation  on  the  assump- 
tion that  the  flutter  condition  is  more  difficult  to  achieve  for  asymmetrical  con- 
• figurations.  Therefore  a procedure  often  used  to  reduce  the  size  of  the  flutter  pres- 

entation task  is  to  analyze  all  configurations  as  being  carried  symmetrically  [ 1 J , £ 2 J , [ 3 J 
i Beyond  this,  in  his  paper  on  "Flutter  of  Aircraft  with  External  Stores”,  presented  at 

the  U.S.  Air  Force  Aircraft/Stores  Compatibility  Symposium  in  19b9,  11. Katz  Indicated 
the  possibility  to  Increase  the  flutter  speed  of  store  configurations  by  "built-in" 
unsymmetries  like  making  the  left-hand  pylon  of  different  stiffness  than  the  right-hand 
pylon.  Of  course  this  procedure  would  only  be  reliable  If  for  a given  aircraft  all 
possible  store  configurations  are  proved  to  flutter  at  higher  speeds  when  the  system 
Is  made  asymmetrical. 


First  indications  of  asymmetrical  store  Instabilities  at  lower  air  speeds  than  meas- 
ured for  the  related  symmetrical  configurations  were  found  during  wind  tunnel  testing  of 
a model,  having  sweepable  wings.  Results  of  this  test  wore  presented  at  the  "Specialists 
Meeting  on  Wing-Wlth-Stores  Flutter"  in  Munich  1974  [4]  . Due  to  the  small  damping  gra- 
dients measured  for  the  asymmetrical  store  configurations, the  determination  of  the  exact 
flutter  speed  was  difficult  and  still  uncertainties  remained  whether  the  measurements 
demonstrate  real  flutter  or  rather  have  to  be  explained  by  forced  vibrations  caused  by 
marginally  stable  modes.  In  the  meantime  a large  number  of  asymmetrical  store  configura- 
tions have  been  investigated  by  analysis  and  further  wind  tunnel  testing  to  establish  the 
physical  background  for  the  flutter  mechanism  of  asymmetrical  stores  and  to  find  out 
whether  unfavourable  effects  of  asymmetries  exist  only  for  mild  flutter  or  could  also 
occur  for  flutter  cases,  exhibiting  large  gradients  of  aerodynamic  damping  with  air- 
speed. Results  of  this  asymmetrical  store  study  will  be  presented  here. 


TECHNICAL  APPROACH 

Representation  of  Aircraft  Structure 

The  aircraft  Investigated  here  features  sweepable  wings  which  require  single 
point  attachments  for  the  Inboard  and  outboard  wing  pylons.  Due  to  this  requirement  the 
values  for  the  attachment  stiffness  in  the  yaw  degree  of  freedom  are  relatively  small 
for  the  wing  and  also  for  the  pylons.  The  flexibilities  of  the  attachments  had  to  be 
considered  by  special  degrees  of  freedom  in  roll,  pitch  and  yaw  for  the  wing  and  in  yaw 
for  the  pylon.  Since  it  was  proved  by  analysis  that  the  flexibility  of  fin,  talleron  and  j 

fuselage  has  no  effect  on  store  flutter, these  components  were  assumed  to  be  rigid  for 
moat  of  the  calculat Ions. 


For  the  representation  of  wing  external  store  dynamics  the  following  generalized 
coordinates  were  Introduced  into  vibration  and  flutter  analyses. 


3 Rigid  aircraft  modes 

4 Elastic  wing  modes  (first  and  second  vertical  bending,  first  torsion,  first 

lateral  bending) 

3 Ming  pivot  modes 

3 Wing  inboard  store  junction  modes 
3 Elastic  Inboard  pylon  modes 
1 Inboard  pylon  pivot  (yaw)  mode 
3 wing  outboard  store  junction  modes 
3 Elastic  outboard  pylon  modes 
1 Outboard  pylon  pivot  (yaw)  mode 


Using  the  Q-R  Algorlthmus,  the  flutter  equation  was  solved  for  the  equivalent  amount 
of  structural  damping  g, necessary  to  provide  harmonical  oscillations.  Structural  damping 
as  measured  in  ground  resonance  tests  has  not  been  considered  here. 


Unsteady  Aerodynamic  Forces 

Three-dimensional  unsteady  aerodynamic  forces  were  calculated  for  the  wing  using  kernel 
function  theories  [5] . According  to  the  operational  speed  of  most  external  store  configu- 
rations, the  calculations  in  general  have  been  performed  for  the  subsonic  region  at  Mach 
number  0.9.  For  comparison  with  flutter  model  test  results  some  calculations  are  based  on 
M * 0.3. 

Investigations  show  that  in  most  cases  the  influence  on  flutter  behaviour  caused  by 
unsteady  aerodynamic  forces  on  external  stores  and  store-wing  interference  effects  is 
small  [6], [7].  Therefore  unsteady  aerodynamic  effects  of  the  store  are  not  considered 
here . 


Mathematical  Description  of  the  Asymmetrical  System 


The  dynamical  behaviour  of  the  symmetrical  clean  aircraft  was  described  by  available 
symmetrical  and  antisymmetrical  component  modes.  Asymmetrical  external  stores  were  added 
to  the  symmetrical  aircraft  by  coupllnq  asymmetrical  store  modes  (representing  the  mass 
and  stiffness  distribution  of  individual  pylon-store  systems  on  each  side  of  the  wing) 
with  the  symmetrical  and  antisymmetrical  clean  aircraft  modes  dynamically.  This  approach 
has  the  advantage  of  using  proved  structural  and  aerodynamical  representations  of  the 
aircraft,  already  established  for  symmetrical  and  antisymmetrical  investigations.  Having 
defined  a set  of  symmetrical  and  antisymmetrical  clean  aircraft  modes  and  asymmetrical  store 
modes, the  vibration  and  flutter  equations  can  be  solved  as  usual. 


Description  of  the  Flutter  Mechanism 


The  mechanism  of  the  symmetrical  and  asymmetrical  store  flutter  on  principle  can  be 
demonstrated  by  Fiq.  1 whtch  shows  the  variation  of  flutter  speed  and  modal  frequencies 
with  increasing  store  radius  of  gyration  for  the  tore  pitch  axis  as  resulted  from  the 
flutter  model  test  mentioned  above. 

The  three  Important  modes  involved  In  the  flutter  mechanism  are: 


First  wing  bending  mode  with  in-phase  lateral  motion  of  the  store 

Store  pitch  mode  exhibiting  large  amplitudes  of  wlnq  torsion 

Store  roll  mode  coupllnq  with  wlnq  bendlnq  and  out-of-phase  lateral  store 
motion 


'LKfW 


TORE  PITCH 


FIG.  1 FLUTTER  SPEED  AND  MODAL  FREQUENCIES  VERSUS  RADIUS  OF  GYRATION 


Frequencies  of  the  store  roll  mode  and  winq  bendinq  mode  are  not  affected  by  variation 
of  the  pitch  Inertia, but  the  distances  to  the  store  pitch  frequency  are  reduced  if  the 
store  is  carried  asymmetrically. 


The  diagram  on  top  Indicates  two  different  flutter  modes  characterized  by  the  coupling 
of  store  pitch  with  wing  bending  for  the  right-hand  branch  and  store  pitch  with  store  roll 
for  the  left-hand  branch. 

Due  to  the  increased  frequency  of  the  wing  bending  mode  for  asymmetrical  store  carriage 
the  value  for  the  pitch  Inertia  related  to  the  minimum  flutter  speed  condition  is  shifted 
to  smaller  values  of  §>,  thus  creating  the  capability  of  asymmetrical  flutter  at  lower  air 
speeds  in  the  range  of  intermediate  values  of  f . 


4 


For  the  atlffneaa  and  maaa  distributions  Investigated  by  thla  diagram, the  pitch-roll 
flutter  of  the  left-hand  branch  seems  to  result  to  higher  flutter  speeds  for  asymmetrical 
stores  Inside  the  range  of  actual  values  ot  f . 

To  confirm  this  flutter  mechanism  as  deduced  from  wind  tunnel  tests  and  to  examine 
the  asymmetrical  flutter  behaviour  under  different  conditions  for  wing  and  pylon  stiff- 
nesses and  store  weights,  an  analytical  study  was  Initiated  which  will  be  discussed  now. 


RESULTS 

Correlation  of  Test  and  Analysis  Results 

For  all  external  store  configurations  Investigated  here,  symmetrical  flutter  was  found 
to  be  more  critical  than  antlsymmetrlcal  flutter.  Therefore  antisymmetr leal  flutter 
will  not  be  considered. 

Fig.  2 shows  a comparison  of  analysis  results  with  test  results  which  have  been  ob- 
tained by  the  flutter  model  with  modified  wing  and  pylon  stiffnesses. 
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FIG.  2 FLUTTER  SPEED  VERSUS  INBOARD  STORE  PITCH  RADIUS  OF  GYRATION  FOR 
SYMMETRICAL  ANO  ASYMMETRICAL  CONFIGURATIONS 
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In  general  the  mechanism  of  flutter  as  found  during  earlier  wind  tunnel  testing 
are  confirmed.  For  symmetrical  store  configurations  the  two  different  branches  of  flutter 
are  separated  well  and  agreement  between  test  and  analysis  Is  sat lsfactory, considering 
that  small  differences  between  structural  model  data  and  design  data  can  effect  the  mini- 
mum flutter  speed  to  occur  at  smaller  values  of  p,  thus  shifting  the  whole  curve  as  indi- 
cated by  the  diagram. 


As  it  turned  out  for  the  configurations  considered  in  this  f lgure,  the  damping  gradi- 
ent In  the  flutter  point  was  very  small  for  asymmetrical  stores  at&Ref.  values  between 
1.2  and  0.8.  This  made  it  extremely  difficult  to  define  the  exact  flutter  speed  by  model 
testing  which  may  explain  the  difference  between  measured  and  calculated  results  for 
asymm  trlcal  stores.  Nevertheless  calculations  and  also  measurements  Indicate  a range 
of  possibly  lower  flutter  speed  for  asymmetrical  store  carriage. 


For  the  left-hand  branch  no  results  are  available  from  test  for  asymmetrical  store 
conf Igurat Ions . Calculated  differences  between  symmetrical  and  asymmetrical  stores  seems 
to  be  less  essential,  but,  as  it  will  be  shown  later,  this  behaviour  could  be  changed  for 
other  combinations  of  wing  and  pylon  stiffnesses  and  store  weights. 


Whether  asymmetrical  store  carriage  will  increase  or  reduce  the  flutter  speed  can 
already  be  deduced  from  changes  of  frequencies  and  mode  shapes  of  corresponding  normal 
modes . 


Fig.  3 demonstrates  the  differences  in  mode  shapes  for  the  modes  important  for  store 
flutter,  depicting  a value  of  store  pitch  inertia  close  to  the  minimum  flutter  speed  con- 
dition for  symmetrical  stores.  Whereas  for  the  symmetrical  configuration  the  pitch  mode 
exhibits  large  wing  torsional  displacements,  being  able  to  generate  severe  flutter, the 
contribution  of  wing  torsion  to  the  asymmetrical  store  pitch  mode  is  small.  The  asymme- 
trical store  pitch  frequency  Is  shifted  below  the  wing  bending  frequency  which  leads  to 
a mild  flutter  case.  The  store  yaw  and  store  roll  modes  In  both  cases  are  well  separated 
from  the  store  pitch  mode  and  are  not  involved  In  the  flutter  mechanism. 


How  the  wing  bending  frequency  and  the  shape  of  the  store  pitch  mode  change  with 
varying  j> -value  Is  shown  by  Fig.  4.  This  fiqure  Illustrates  for  the  symmetrical  configu- 
ration the  tuning  of  frequency  between  store  pitch  and  wing  bending  and  the  variation  of 
the  pitch  nodal  line  from  the  most  critical  position  ati^pRef.  ■ 1.2  to  a less  critical 
forward  position  atS)^  rc{  * 0.8.  For  the  asymmetrical  configuration  the  worst  nodal  line 
position  is  obtained  for  a$^RCf.  value  of  1.0, whereas  for  the  symmetrical  store  carriage 
the  nodal  line  Is  already  shifted  forward  for  this  condition. 


Following  the  experience  gained  from  analyses  and  model  testing, dif ferenees  in  the 
flutter  behaviour  of  symmetrical  and  asymmetrical  configurations  can  be  affected  by 

. reduction  of  generalized  mass  for  the  asymmetrical  case  which  is  easier 

to  excite  and  therefore  may  lower  the  flutter  speed 

. change  of  mode  shapes  on  the  unloaded  wing  side,  which  at  frequencies  of  the 

three  store  modes  shows  contributions  of  the  wing  bending  mode  only.  The  loaded 

wing  side  only  exhibits  wing  torsional  motions  in  the  pitch  mode,  capable  to 
initiate  flutter.  Energy  transfer  from  the  loaded  side  to  the  unloaded  side 
therefore  results  to  a more  stable  overall  configuration 

. changed  nodal  line  position  of  the  pitch  mode  which  seems  to  be  the  most  im- 
portant effect  of  unsymmetries  because  this  can  shift  the  minimum  flutter 
speed  to  a condition,  at  which  the  symmetrical  configuration  may  be  less 
critical 

. Increased  wing  bending  frequency  and  decreased  store  roll  frequency  for  the 
asymmetrical  case  which  can  create  a coupling  with  the  store  pitch  mode  by 
smaller  aerodynamic  forces  than  necessary  for  symmetrical  configurations  to 
obtain  the  flutter  condition. 

The  differences  between  symmetrical  and  asymmetrical  store  flutter  seems  to  be  re- 
duced for  mild  flutter  cases  if  actual  values  of  structural  damping  as  measured  in  ground 
resonance  tests  are  considered. 


Analytical  Results  for  Different  Stores  Carried  on  Inboard  and  Outboard  Wing  Pylons 

The  flutter  behaviour  will  now  be  discussed  in  detail  for  different  inboard  and  out- 
board store  conditions. 
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One  of  the  largest  stores  usually  being  carried  on  the  wing  Is  the  external  fuel  tank 
For  most  aircrafts  the  lowest  flutter  speed  will  be  obtained  for  the  maximum  fuel  con- 
dition, representing  large  values  of  store  weight  and  store  radius  of  gyration  which  gen- 
erate flutter  according  to  the  "right-hand"  flutter  branch. 


FIG.  5 ASYMMETRICAL  FUEL  EMPTYING  OF  INBOARD  PYLON  TANKS 


LEFT  HAND  TANK  100%  and  50%  FUEL 
RIGHT  HAND  TANK  0 100%  FUEL 


M - 0 0.  WET  WING.  A. > 2S  , PYLON  STIFFNESS  00%  OF  NOMINAL 


Fig.  5 demonstrates  the  variation  of  flutter  speed  with  asymmetrical  fuel  emptying 
of  the  inboard  pylon  tanks  at  a wing  sweep  position  of  25  degrees.  The  internal  wing 
tanks  are  considered  to  be  full  because  in  normal  flight  operation  the  external  tanks 
are  emptied  before  the  wing  fuel  is  being  consumed.  For  the  lower  curve  one  external 
tank  is  assumed  to  contain  maximum  fuel  whereas  the  fuel  of  the  other  tank  varies  from 
zero  to  lOO  ».  The  lowest  flutter  speed  is  reached  for  the  fuel  state  with  both  tanks 
being  full  which  matches  about  the  minimum  flutter  speed  condition  for  this  configuration 
Emptying  only  one  of  the  tanks,  the  wing  bending  frequency  lies  above  the  store  pitch 
frequency  which  provides  a very  mild  flutter  at  considerably  higher  flutter  speeds. 


The  upper  curve  demonstrates  the  flutter  Trend  of  the  tanks  when  the  fuel  state  of 
one  tank  is  kept  in  the  50  % condition  and  the  other  is  being  varied  from  zero  to 
lOO  *.  Two  different  flutter  modes  are  obtained,  a very  mild  flutter  at  4.1  Hz  when  at 
least  one  tank  is  about  full  and  a more  excessive  flutter  at  6.2  Hz  which  is  dominated 
by  the  pitch  mode  of  the  tank  with  50  * fuel  content. 

For  this  configuration  it  is  proved  by  the  results  shown  in  the  diagram  that  all  symmet- 
rical and  asymmetrical  fuel  state  conditions, which  may  occur  during  operational  flight 
or  may  result  from  fuel  emptying  failure  cases, are  covered  by  the  flutter  speed  of  the 
symmetrical  configuration  with  full  tanks. 


How  the  flutter  behaviour  changes  if  outboard  stores  are  added  to  these  tank  configu- 
rations is  illustrated  by  the  next  figures. 


*\  A *s  \S  A«. 

t V 


KH»M1  M H MU  i\’NI|  N! 


FIG  6 


FIG  1 


ItH  HANOI  ANA  100%  Kill 
HK.H1  HANOI  ANA  0 100% 


UM  HANOI  ANA  H0%»Utt 
KIliHI  MANO  1 ANA  0 100%  HUt 


M 0«  Mil  MINU.  fVl  ON  SI  It  HlltM  *0%  Of  NOMINAL 


m oo  mnioNvi,  r*  ion  *io»niss»o%o»  nominal 


'"1 


ASYMMETRICAL  FUEL  EMPTYING  OF  INBOARO  PYLON  TANKS  WITH  SYMME  TRICAl 
OUTBOARD  STORES 


In  FI>|.  *>  the  leant  b am  piesented  which  have  been  achieved  foi  two  dltteient  wing 
■weep  positions.  One  external  tank  Ih  assumed  to  h«>  full  and  the  otlvei  tank  va>  le« 
from  ien>  to  too  » fuel  state.  Fot  t tie  2k  degree  wing  sweep  post!  ton  t lie  highest  fluttei 
speed  was  oaloulat e»l  for  ttie  symmetrical  full/full  condition.  Kmptylng  one  tank  only  t he 
fluttei  speed  is  decrees  Ing  amt  reaches  ttie  lowest  value  at  t tie  empty  full  condition. 

This  behaviour  Is  not  changed  If  a relatively  large  value  of  2.k  t Htiuctuial  damptnu  Is 
considered, which  proves, that  the  damping  gradient  In  the  tluttei  point  Is  compaiable  foi 
symmetilcal  and  asymmetrical  fuel  condittona. 

This  trend  changes  to  the  opposite  If  the  wing  Is  swept  hack  . In  tills  case  ttie  symmet 
rlcal  full/full  condition  represents  the  woiBt  condition  and  would  deal  also  ttie  tluttei 
of  asymmetrical  fuel  states.  The  second  fluttei  mode  at  h.O  tta  Is  not  tmpoitant  tot  tluttei 
because  It  Indicates  hlghei  flutter  speeds  for  all  fuel  states. 


Fig.  7 shows  the  flutter  speeds  as  resulted  foi  the  same  store  oonf lgut at  ton  it  one 
tank  Is  being  kept  In  the  lO  » fuel  condition.  For  the  2k  degiee  wing  sweep  position  ttie 
results  ale  slmllai  to  thosi  obtained  foi  the  tank  oonf tgui at  Ion  without  outboard  atotea 
(see  Fig.  k)  . t'onsldeilng  2.1  * Htiuctuial  damping,  t lie  lowest  fluttei  speed  occm  s foi 
the  symmetrical  condition. 

At  <k  degree  wing  sweep  the  lowest  fluttei  speed  Is  achieved  foi  the  ko  t tin'  * fuel 
condition  which  was  alteady  Included  In  the  diagram  shown  hefoie. 

It  Is  striking  that  all  lesults  obtained  foi  ttie  tank  fuel  emptying  reveal  only  little 
effect  of  asymmetries  on  the  second  fluttei  mode  which  Is  dominated  by  ttie  stole  pitch 
mode  of  the  tank  with  ko  * fuel,  coupling  with  ttie  wing  bending  mode  at  a frequency  ot 
about  h IU.  Tills  Is  also  evident  from  ttie  fluttei  speed/damptng  curves  foi  ttie  tuboaid 
store  carriage  shown  In  Fig.  8. 


FIG  8 FIG  9 

DAMPING  OF  CRITICAL  FLUTTER  MODES  FOR  DIFFERENT 
INBOARD  TANK  FUEL  CONDITIONS 
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The  b.2  tlz  flutter  movies  are  closely  spaced  for  all  possible  tank  combinations  where- 
as the  4.1  Hz  flutter  movies  exhibit  a wlvie  scatter  of  the  flutter  speed  and  also  of  the 
damping  gradient.  This  phenomenon  Indicates  that  the  effect  of  asymmetries  on  flutter  is 
large  if  the  considered  store  radius  of  gyration  Is  close  to  the  value  which  defines  the 
minimum  flutter  speed,  figure  8 also  demonst rates  for  this  special  store  configuration 
that  in  agreement  with  common  flutter  philosophy  the  flutter  speed  calculated  for  the 
symmetrical  configuration  would  be  safe  for  all  possible  asymmet r leal  conditions,  not 
dependent  on  structural  dampings  measured  for  this  aircraft  in  ground  resonance  test. 

But  deviation  from  this  rule  is  proved  by  Fig.  4 which  shows  the  or  1 1 lea  1 flutter  movies 
for  the  same  configurations  but  carrying  in  addition  stores  on  the  outboard  wing  stations. 
In  this  case  the  flutter  speed  resulted  for  the  full/full  condition  would  not  be  safe  for 
the  asymmetrical  empty/full  tank  configuration  which  flutters  at  considerably  lower  air- 
speed, indicating  about  the  same  decrease  of  damping. 


The  favourable  or  unfavourable  effect  of  asymmet ry,  dependent  on  the  fact  whether 
additional  stores  are  carried  on  the  outboard  wing  station,  can  be  explained  by  the  de- 
crease of  the  wing  bending  frequency  due  to  the  outboard  store.  For  the  symmetrical  full 
tank  configuration  this  effect  is  beneficial  because  it  detunes  the  wing  bending  and 
store  pitch  frequencies  whereas  for  the  asymmetrical  tank  this  frequencies  are  tuned  by 
the  reduction  of  the  wing  bending  frequency.  Now  also  for  asymmetrical  configurations  the 
wing  bending  frequency  is  below  the  store  pitch  frequency  which  creates  the  capability 
of  severe  flutter. 


The  fol lowing  example  Illustrates 
attached  by  adapter  which  reduces  the 
stiffness  the  cons Idered  store  radius 
expected  to  match  the  minimum  flutter 
when  the  wing  is  swept  to  the  forward 


the  asymmetrical  flutter  of  an  inboard  wing  store 
stiffness  of  the  pylon.  Due  to  this  reduction  in 
of  gyration  for  pitch  of£/jpR.,{t  » 0.47  can  be 
speed  condition  for  asymmetrical  store  carriage 
pos i t ion . 
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VARIATION  OF  FLUTTER  SPEED  WITH  WING  SWEEP  ANGLE  FOR  SYMMETRICAL  AND  ASYMMETRICAL 
INBOARD  STORE  CARRIAGE  WITH  AND  WITHOUT  SYMMETRICAL  OUTBOARD  STORES 

M * 01.  DRY  WING  PVlON  STIFFS!  S$  10\  OF  NORMNAi 


The  expected  behaviour  in  confirmed  by  the  results  for  the  inboard  store  carrlaqe, 
depicted  in  Fig.  10  for  different  wing  sweep  positions.  For  25  degrees  wing  sweep  flutter 
occurs  at  considerably  lower  flutter  speed  for  the  asymmetrical  configuration.  This 
diagram  also  demonstrates  the  largo  influence  of  the  wing  sweep  angle  on  the  asymme- 
trical flutter.  At  wing  sweep  positions  higher  than  35  degrees  asymmetrical  flutter  is 
not  Important. 


The  flutter  behaviour  is  changed  completely  if  additional  stores  are  mounted  on  out- 
board wing  pylons.  Fig.  11  indicates  very  mild  flutter  at  small  and  intermediate  sweep 
angles  for  the  symmetrical  and  also  for  the  asymmetrical  store  configuration.  The  criti- 
cal flutter  mode  at  large  wing  sweep  angles  is  dominated  by  the  outboard  store  modes  and 
is  therefore  not  effected  by  asymmetries  of  the  inboard  stores. 


The  flutter  critical  situation  for  the  Inboard  store  configuration  at  25  degree  sweep 
angle  is  depicted  in  Fig.  12.  The  shapes  of  the  wing  bending  and  store  pitch  mode,  which 
are  important  for  this  flutter  case, arc  about  identical  for  the  wing  which  carries  the 
store.  But  frequency  separation  is  better  for  the  symmetrical  configuration.  From  the 
damping  and  frequency  plot  of  the  four  flutter  modes  up  to  10  Hr  it  can  be  seen  that 
flutter  occurs  for  both  configurations  by  the  coupling  between  wing  bending  and  store 
pitch.  Due  to  the  higher  wing  bending  frequency  at  zero  air  speed  the  coupling  point  is 
reached  at  lower  flutter  speed  for  the  asymmetrical  store.  It  is  also  evident,  that  store 
yaw  (mode  No.  2)  and  store  roll  (mode  No.  4)  are  not  involved  in  the  flutter  mechanism. 


The  next  flutter  case  which  will  be  discussed  now  is  a more  complicated  example  for 
the  "left-hand"  flutter  branch,  characterised  by  the  store  pitch/store  roll  flutter. 

Fig.  13  shows  the  four  normal  modes  which  could  contribute  to  the  store  flutter.  For 
both  configurations,  the  symmetrical  and  also  the  asymmetrical,  the  mode  shapes  are 
marked  by  a strong  coupling  between  store  yaw,  pitch  and  roll,  due  to  closely  spaced 
frequencies  of  the  store  component  modes.  The  second  and  also  the  third  mode  indicate 
contributions  of  store  pitch, but  the  second  mode  show  much  larger  motion  of  store  roll/ 
wing  bending, wheroas  the  twist  nodal  line  of  the  wing  exhibits  a less  critical  forward 
position. 
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It  can  be  assumed,  therefore,  that  the  second  mode  will  behave  more  like  a store 
roll/wlnq  bending  mode  and  the  third  mode  will  be  the  pitch  mode  which  becomes  critical. 
Differences  In  mode  shapes  and  frequencies  between  the  symmetrical  and  the  asymmetrical 
configuration  are  very  small.  In  this  case  the  effect  of  energy  transfer  between  both 
sides  of  the  wing  can  be  expected  to  result  to  a higher  flutter  speed  for  the  asymmetrical 
condition. 


SYMM.  STORE 


-777-  • 1.08  -4~. 


ASYMM.  STORE 


0,71 


FIG.  14  DAMPING  AND  FREQUENCY  VERSUS  FLUTTER  SPEED  FOR  SYMMETRICAL  AND 
ASYMMETRICAL  INBOARD  STORE  CARRIAGE 
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The  damping  and  frequency  plot  of  Fig.  14  prove  that  these  assumptions  are  true.  For 
both  configurations  the  coupling  between  the  second  and  third  flutter  mode  defines  the 
lowest  flutter  point.  Due  to  the  fact  that  the  store  roll  frequency  Is  below  the  store 
pitch  frequency,  the  flutter  behaviour  of  the  symmetrical  configuration  is  rather  exces- 
sive. As  expected,  the  asymmetrical  store  flutters  at  much  higher  air  speeds,  because 
neither  the  frequencies  nor  the  wing  twist  nodal  line  are  changed  to  a worse  condition 
by  the  asymmetry. 


The  next  configuration  considered  deals  with  a single  pylon  store  on  outboard  wing. 
It  is  quite  clear,  that  flutter  trends  established  for  Inboard  stores  to  Indicate  the 
effect  of  asymmetries  can  not  be  used  directly  for  outboard  stores.  It  has  to  be  con- 
sidered that  the  Influence  of  outboard  stores  on  the  wing  bendlnq  frequency  is  much 
larger  but  in  general  similar  effects  of  asymmetries  can  be  expected. 


The  first  two  modes  of  Fig.  15,  showing  the  normal  modus  for  this  configuration, 
indicate  a strong  coupling  between  wing  bending  and  store  yaw  for  the  asymmetrical  case 
due  to  the  increased  frequency  of  wing  bending.  For  both,  the  symmetrical  and  the  asym- 
metrical configuration,  two  modes  with  large  store  pitch  motions  (modes  No.  3 and  4)  are 

evident  which  are  generated  by  a coupling  between  store  pitch  and  wing  yaw.  Both  modes 

are  able  to  initiate  flutter  but  the  mode  with  the  more  critical  lower  frequency  exhibits 

also  the  more  critical  wing  nodal  line  position.  Because  of  the  smaller  distance  between 

the  wing  bending  and  the  store  pitch  frequency  the  asymmetrical  store  must  be  expected 
to  flutter  at  the  lower  air  speed.  Due  to  the  very  high  frequency  (about  14  Hz)  the  store 
roll  mode  is  not  involved  in  the  flutter  mechanism  and  is  therefore  not  shown  here. 
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FIG.  16  DAMPING  AND  FREQUENCY  VERSUS  FLUTTER  SPEED  FOR  SYMMETRICAL  AND 
ASYMMETRICAL  OUTBOARD  STORE  CARRIAGE 

M * 0 9 DRY  WING. A « 45  PYLON  STIFFNESS  70%  OF  NOMINAL 

The  variation  of  dampings  and  frequencies  are  demonstrated  for  this  outboard  store 
configuration  in  Fig.  16.  The  specific  behaviour  of  this  light  weight  store  is  reflected 
by  the  flutter  modes  No.  1 and  2 of  the  symmetrical  configuration  which  show  that  these 
modes  are  interchanging.  Finally  the  second  flutter  mode  is  dominated  by  wing  bendinq 
and  the  frequency  crossing  of  this  mode  and  the  store  pitch  mode  defines  the  flutter  point 

For  the  asymmetrical  store  small  aerodynamic  forces  are  sufficient  to  sepaiate  wing 
bending  and  store  yaw  (Modes  No.  1 and  2).  In  this  case  the  wing  motion  of  the  store 
pitch  mode  (Mode  No.  3)  generates  much  smaller  aerodynamic  forces,  but  due  to  the  higher 
wing  bending  frequency  the  critical  flutter  condition  is  reached  at  lower  flutter  speed 
than  obtained  for  the  symmetrical  configuration. 

Not  considering  the  more  complicated  coupling  between  the  modes  engaged,  this  example 
shows  that  the  fundamental  mechanism  of  asymmetrical  store  flutter  is  similar  for  in- 
board and  outboard  stores. 
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Effect  of  Asymmetrical  Stlf  fne  s a D_1  at  r lbut  Ion 

The  last  figures  presented  here  demonstrate  the  effect  of  asymmet r leal  pylon  stiff- 
nesses on  the  flutter  speed.  Two  different  Inboard  store  conditions  are  investigated. 

For  the  first  store  a radius  of  gyration  was  assumed  which  matches  the  minimum  flutter 
speed  condition  for  symmetrical  store  carriage  and  which  results  to  higher  flutter  speeds 
if  only  one  store  is  being  carried  (Pig.  17(,  for  the  second  store,  attached  to  the 
pylon  by  an  adapter,  a -value  was  chosen  which  is  close  to  the  minimum  flutter  speed 
condition  if  the  store  is  carried  asymmetrically  (Fig.  181.  For  this  case  the  more 
critical  flutter  was  found  for  the  asymmetrical  configuration  with  one  store  only. 


RIGHT  SIDE  PI  ION  STIFFNESS  (•/.) 


FIG. 17 


FIG. 18 


VARIATION  OF  FLUTTER  SPEED  WITH  ASYMMETRY 
OF  INBOARD  PYLON  STIFFNESS 


M.OI  ORV  WING  A- T5MEfT  NANO  AVION  STIFFNtSS  >00% 


In  Fig.  17  and  18  the  flutter  speeds  are  plotted  versus  stiffness  of  cv>-  pylon  where- 
as the  stiffness  of  the  other  pylon  is  being  kept  at  lOO  1.  Both  figures  - llcate  a change 
of  the  critical  flutter  mode  when  the  stiffness  of  one  pylon  decreases  of  the 

nominal  stiffness. 

When  the  stiffness  of  one  pylon  is  increased  to  1201,  a further  flutter  mo  ar  found 
at  very  high  flutter  speed  outside  the  range  of  the  diagrams.  From  this  the  ..elusion 
can  be  drawn  that  each  wing  side  creates  its  corresponding  flutter  case,  each  navlng  a 
minimum  flutter  speed  at  defined  values  of  pylon  stiffness.  For  the  6.0  Hz  cur  c in 
Fig.  17,  which  refers  to  the  flutter  of  the  pylon  with  lOO  » stiffness,  the  minimum 
flutter  speed  is  reached  at  10O  » stiffness,  corresponding  to  the  f-value  which  already 
matches  the  minimum  flutter  speed  condition  at  nominal  stiffness.  For  the  6.6  Hz  flutter 
mode  of  Fig.  18  the  minimum  is  not  reached  yet,  which  is  reflected  by  the  decreasing 
flutter  speed  when  the  stiffness  of  emt  pylon  is  being  Increased.  Both  diagrams  also 
indicate  that  for  the  flutter  mode  related  to  smaller  stiffnesses  the  flutter  speed  is 
still  reducing  for  further  reduction  of  the  stiffness. 


The  results  of  this  calculations  demonstrate  that  asymmetry  of  pylon  stiffness  can 
be  very  important  for  values  of  store  inertia  and  pylon  stiffness,  which  are  close  to 
the  minimum  flutter  speed  condition.  For  the  cases  considered  here  the  Influence  is 
reduced,  if  structural  damplnq  is  taken  into  account.  For  final  judgement  this  phenome- 
non needs  to  be  further  investigated,  especially  in  view  of  stores  having  intermediate 
and  smaller  radii  of  gyration. 

Differences  in  pylon  stiffnesses  inside  the  usual  range  of  small  tolerances  seems  to 
be  less  important  for  flutter  but  can  create  considerable  problems  during  the  evaluation 
of  flight  test  data,  caused  by  closely  spaced  frequencies,  which  may  result  to  misleading 
high  damping  values.  In  this  case  a verification  of  test  results  by  flutter  calculations 
is  even  more  Important  and  it  may  be  necessary  to  provide  the  flutter  clearance  of  store 
configurations  with  marginal  flutter  speeds  rather  by  analysis  than  by  flight  test. 


CONCLUSION 

Coming  to  the  conclusion  it  should  be  kept  in  mind  that  all  results,  demonstrated 
here,  have  been  obtained  for  wing  pylons  with  single  point  attachments,  due  to  the  require- 
ment of  sweepable  wings.  For  other  pylon  attachments,  changes  of  the  critical  flutter 
conditions  by  asymmetrical  store  carriage  may  be  different.  Nevertheless,  some  general 
rules  can  be  established  by  the  results  of  this  study  which  indicate  the  trend  of  flutter 
changes  caused  by  asymmetries. 

The  fundamental  mechanism  of  store  flutter  is  characterized  by  the  aerodynamic  coupling 
of  the  store  pitch  mode,  producing  large  wing  torsional  motions,  with  wing  bending.  Two 
modes  with  large  wing  bending  motions  are  existent  which  are  defined  by  the  "in-phase“ 
and  *out-of-phase“  coupling  with  lateral  store  motions.  According  to  this,  two  different 
flutter  cases  had  to  be  considered  which  are  able  to  generate  low  flutter  speeds  at  very 
large  or  very  small  values  of  store  inertia  about  the  pitch  axis.  For  intermediate  values 
of  store  inertia  the  influence  of  asymmetries  is  less  important  because  reasonable  high 
flutter  speeds  can  be  expected  for  this  region. 

It  has  been  shown  by  this  study  that  asymmetrical  store  carriage  or  asymmetrical 
stiffness  distribution  results  to  changes  of  the  wing  bending  frequency  and  the  wing 
nodal  line  position  of  the  store  pitch  mode,  which  can  generate  either  lower  or  higher 
flutter  speeds.  For  store  configurations  with  values  of  pylon  stiffness  and  store  inertia 
close  to  the  minimum  flutter  speed  conditio1'  the  changes  by  asymmetries  can  be  most  effec- 
tive, caused  by  the  tuning  or  detuning  of  wing  bending  and  store  pitch  frequencies.  It 
also  has  been  found  that  the  influence  of  asymmetries  on  flutter  Is  beneficial,  if  the 
change  in  nodal  line  position  or  wing  bending  frequency  is  small.  In  this  case  energy 
transfer  between  the  more  critical  and  the  less  critical  wing  is  decisive. 

Knowing  the  flutter  mechanism  it  is  quite  clear  that  changes  of  flutter  parameters, 
like  store  weight,  wing  sweep  position  and  the  carriage  of  additional  stores  on  he  out- 
board wing  pylon  which  influence  the  wing  bending  frequency,  or  changes  of  pyloi  stiff- 
nesses which  effect  the  store  pitch  frequency,  will  also  change  the  effect  of  asymmetries 
considerably  as  it  was  shown  by  the  results. 

The  usual  procedure  to  analyze  all  store  configurations  as  being  carried  symmetrically 
can  not  be  followed  for  all  aircrafts.  This  does  not  mean  that  all  possible  asymmetrical 
configurations  have  to  be  investigated. 

It  is  recommended  to  establish  flutter  trends  by  variation  of  important  parameters 
before  actual  store  configurations  are  being  calculated.  Once  the  regions  with  possibly 
lower  flutter  speeds  of  asymmetrical  stores  are  def ined,  those  configurations  can  be 
selected  which  have  to  be  investigated. 
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SUMMARY 

A wind  tunnel  test  program  scheduled  for  completion  in  April  1978  is  being  sponsored  by  the  Air  Force 
Flight  Dynamics  Laboratory  for  demonstration  of  active  wing/store  flutter  suppression  systems  on  a light- 
weight fighter  aircraft.  Northrop  Corporation  was  selected  as  contractor  of  the  program  which  included 
preliminary  design,  final  design,  fabrication  and  testing  of  a wind  tunnel  model.  The  present  paper 
presents  preliminary  results  of  the  design  analysis  and  the  test  program  which  was  conducted  at  the  NASA/ 
Langley  Transonic  Dynamics  Tunnel.  Three  configurations  were  selected  for  final  testing.  Two  of  these 
configurations  were  deliberately  designed  to  exhibit  low  flutter  speeds  with  rapid  reductions  in  damping 
at  the  incipient  flutter  condition.  After  initial  tunnel  entries  which  showed  the  need  for  certain  im- 
provements in  the  model  and  the  control  system  design,  substantial  increases  in  the  flutter  speeds  were 
achieved  using  both  leading  and  trailing  edge  control  surfaces  separately.  For  the  most  critical  config- 
uration a demonstrated  improvement  of  18%  and  a projected  improvement  of  29%  in  the  dynamic  pressure  were 
accompl ished . 


1.  INTRODUCTION 

The  increase  in  strength  of  materials  and  the  use  of  thin,  low-drag  airfoils  have  led  to  flexible 
and  more  flutter-prone  wings.  With  many  combinations  of  external  stores,  modern  fighter  aircraft  have  a 
wide  variety  of  flutter  placards  that  restrict  operational  use.  In  fact,  an  aircraft  carrying  external 
stores  may  have  several  flutter  speed  restrictions  on  a single  flight.  In  order  to  reduce  the  severity 
of  these  restrictions  and  take  advantage  of  recent  improvements  in  controls  technology  for  high-performance 
aircraft,  a logical  approach  is  to  develop  a flutter  suppression  system  that  eventually  could  employ  the 
same  components  used  for  conventional  stability  augmentation.  Before  a working  flutter  suppression  sys- 
tem can  be  installed  in  an  aircraft,  a number  of  design  aspects  have  to  be  considered.  These  aspects 
include  selection  of  control  surfaces  and  related  actuation  systems,  definition  of  appropriate  control 
laws  applicable  to  a large  variety  of  configurations  and  flight  conditions,  and  development  of  redun- 
dancy and  fail-safe  features. 

Automatic  feedback  systems  controlling  the  aeroelastic  response  have  been  used  successfully  on 
several  aircraft.  The  applications  to  date  have  been  for  augmentation  of  static  stability,  improvement 
of  ride  quality,  reduction  of  maneuver  loads  and  suppression  of  structural  loads  induced  by  atmospheric 
turbulence.  Considerable  interest  has  emerged  in  recent  years  towards  adding  an  active  flutter  control 
function  to  the  flieht  control  systems.  A pioneering  investigation  was  performed  by  the  Boeing  Company 
and  Honeywell,  Inc.'  paving  the  way  for  a flight  test  program  that  employed  an  Air  Force  Flight  Re- 
search B-52  to  demonstrate  the  feasibility  of  active  flutter  suppression.  In  addition,  active  flutter 
suppression  tests  were  conducted  on  a wind  tunnel  model  of  the  B-52  to  obtain  comparisons  with  flight 

test  data'  Other  efforts  within  the  United  States  included  a wind  tunnel  test  program^)  on  an  SST-  ; 

type  wing  and  an  analytical  study^*5)  of  a wing/store  flutter  suppression  system.  Outside  of  the  United 
States,  considerable  activity  in  the  area  of  flutter  suppression  has  been  evidenced  by  several  publica- 
tions from  England,  France  and  Germany.  Significant  contributions  in  terms  of  wind  tunnel  test  demon- 
strations include  References  (6)  and  (7).  An  interesting  analytical  study  of  active  flutter  suppression  is 
presented  in  Reference  8. 

In  March  1976  the  Air  Force  Flight  Dynamics  Laboratory  (AFFDL)  initiated  the  sponsorship  of  a wind 
tunnel  test  demonstration  of  active  wing/store  flutter  suppression  systfems  for  a lightweight  fighter  air- 
craft. Northrop  Corporation  was  awarded  the  contract  which  has  been  accomplished  in  four  phases:  pre- 
liminary design,  final  design,  fabrication  and  testing  of  a wind  tunnel  model.  The  present  paper  pre- 
sents the  design  analysis  and  the  test  program,  which  included  three  entries  at  the  NASA/Langley  Trans- 
onic Dynamics  Tunnel. 

2.  DESCRIPTION  OF  THE  MODEL 

The  test  specimen  is  a 30%  scale,  semi -span  (right  half)  model  of  a lightweight  fighter  aircraft. 

The  configuration  is  characterized  by  a moderately  swept  wing  with  a large,  highly  swept  leading  edge 
extension  (LEX)  at  the  root,  differential  area  ruling  of  the  fuselage  and  under-wing  engine  inlets  with 
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slots  for  fuselage  boundary  layer  diversion.  The  model  was  designed  such  that,  for  three  selected  store 
configurations,  the  unaugmented  flutter  speeds  plus  the  desired  Improvements  with  the  AFSS  (Active  Flutter 
Suppression  Systems)  operative,  could  be  demonstrated  within  the  test  limits  of  the  NASA/Langley  16-foot 
Transonic  Dynamics  Tunnel  (TDT) . Pertinent  scale  factors  based  on  freon  as  the  flow  medium  are  presented 
In  Table  1. 

The  model  simulates  all  Important  wing,  pylon,  rigid  body,  and  fuselage  degrees  of  freedom  that  are 
required  to  provide  correct  modal  coupling  for  flutter.  A schematic  drawing  is  shown  In  figure  1.  The 
wing  has  an  aluminum  spar  with  segmented  balsa  wood  covering.  The  half  fuselage  Is  simulated  by  a seg- 
mented metal  shell  attached  to  a magnesium  beam  which  Is  restrained  In  Its  lateral  degrees  of  freedom  by 
three  bar  mechanisms  mounted  on  the  tunnel  wall.  In  addition,  the  model  Is  supported  by  a lift  cable  and 
by  two  preloaded  (450N)  fuselage  cables  running  around  pulleys  so  that  the  model  Is  free  to  move  In 
plunge  and  pitch  and,  to  a lesser  degree.  In  the  axial  direction.  To  obtain  smooth  flow  past  the  fuselage, 
a splitter  plate  Is  Installed  close  to  the  symmetry  plane.  The  model  is  trimmed  by  an  all-movable  and 
remotely  controlled  horizontal  tall.  A narrow-span  leading  edge  flap  and  a trailing  edge  surface  were 
selected  as  the  active  flutter  suppression  devices.  They  are  actuated  by  miniature  hydraulic  actuators 
designed  and  supplied  by  the  Boeing  Company,  Wichita  Division.  Accelerometers  are  Installed  Inside  the 
wing  to  sense  the  vibration.  After  proper  conditioning,  the  signals  are  fed  to  the  servo-valve  which 
supplies  hydraulic  pulses  to  the  actuator.  A potentiometer  is  installed  in  the  hinge  mechanism  to  sense 
the  control  surface  rotation  and  provide  Input  to  the  actuator  servo. 


DIMENSION 

SYMBOL 

SCALE 

Length 

L 

0.300 

Velocity 

L/T 

0.451 

Time 

T 

0.665 

Frequency 

l/T 

1.503 

Dynamic  Pressure 

M/LT2 

0.127 

Mass 

M 

0.0169 

Force 

ML/T2 

0.0115 

Flexibility 

t2/m 

26.18 

Table  1.  Model  Scale  Factors 

3.  ANALYTICAL  STUDY 

3.1  Initial  Flutter  Analysis 

Several  flutter  analyses  were  made  to  determine  the  most  interesting  configurations  for  the  test 
program.  Trend  studies  were  performed  with  variations  of  store  mass,  center  of  gravity,  moment  of  in- 
ertia, pylon  attachment,  etc.  Linearized  finite  element  models  of  the  aeroelastlc  system  were  employed 
In  the  conventional  manner  with  the  dynamic  equation  formulated  in  terms  of  generalized  coordinates  in- 
cluding a representative  number  of  natural  vibration  modes  and  appropriate  rigid  body  modes.  Unsteady 
aerodynamic  loadings  were  computed  by  use  of  the  doublet-lattice  method  at  three  discrete  Mach  numbers 
In  the  subsonic  regime.  Three  external  store  configurations  that  differ  substantially  in  the  flutter 
characteristics  were  selected  for  further  analysis  and  wind  tunnel  testing.  The  configurations  are: 

Configuration  (A):  Tip  launcher  rail:  AIM-9E  (Sidewinder) 

Tip  pylon  (95Z  span):  Not  Installed 

Inboard  pylon  (65?  span):  AIM-7  Sparrow  (3"  aft) 

Configuration  (B) : Tip  launcher  rail:  Empty 
Tip  pylon:  AIM-7  (3"  aft) 

Inboard  pylon:  Not  installed 

Configuration  (C) : Tip  launcher  rail:  Empty 

Tip  pylon:  AIM-9E  (6”  aft) 

Inboard  pylon:  Not  installed 

Figure  2 presents  analytical  results  for  three  selected  configurations.  The  preliminary  flutter 
analysis  of  configuration  (A)  predicted  a "hump"  mode  at  about  13  Hz  in  model  scale  with  a slow  decrease 
In  danqplng  as  the  flutter  speed  Is  approached.  As  shown  in  Figure  2,  this  characteristic  results  In  an 
unusual  Increase  of  the  flutter  dynamic  pressure  at  Mach  numbers  approaching  unity.  Configuration  (B) 
exhibits  a conventional  bending-torsion  type  of  flutter  at  around  7 Hz  with  a violent  onset.  It  should 
be  noted  that  this  flutter  condition  was  created  by  attaching  a wing  pylon  in  the  tip  region  with  the 
store  located  In  an  extreme,  aft  position.  Similarly,  a violent  type  of  flutter  was  predicted  also  for 
Configuration  (C),  but  (due  to  the  reduced  store  mass)  at  the  higher  frequency  10.5  Hz. 

3.2  Design  Analysis 

Design  analyses  of  the  active  flutter  suppression  systems  for  each  configuration  were  performed  to 
establish  control  laws,  compressibility  effects,  transducer  and  control  surface  locations,  etc.  For  this 
purpose  new  analytic  techniques  were  developed  at  Northrop  to  cope  with  the  frequency-dependence  of  the 
control  syster  terms,  which  render  the  conventional  eigenvalue  solution  techniques  impractical.  One 


satisfactory  solution  method  used  In  the  present  program  Is  the  Characteristic  Diagram  technique  sug- 
gested by  Landahl^  and  developed  by  Northrop  Ov)  jnto  an  operational  program.  The  method  Involves 
plotting  of  the  characteristic  equation  in  the  complex  plane  at  a given  speed/altitude  combination  with 
the  frequency  as  the  parameter.  Input  frequency  Is  specified  by  a subroutine  that  computes  an  automatic 
frequency  sweep.  The  sweep  rate  is  reduced  as  a root  Is  approached,  thus  defining  the  characteristic 
diagram  with  sufficient  accuracy.  The  damping  of  the  dynamic  system  is  determined  directly  from  the 
shape  of  the  diagram. 

Another  analysis  technique,  also  developed  at  Northrop,  is  the  Transfer  Function  Synthesis  (TFS) 
technique^11^  which  has  proven  to  be  both  practical  and  accurate  for  most  of  the  present  design  analyses. 
The  TFS  technique  is  a computerized  procedure  for  obtaining  closed-form  approximations  to  aeroelastlc 
transfer  functions.  Including  the  effects  of  inertial,  elastic  and  aerodynamic  forces.  The  transfer 
functions  are  synthesized  from  analytical  frequency  response  data  as  ratios  of  polynomials  in  the  La- 
place variable  S.  Figure  3 presents  typical  computed  vertical  and  angular  acceleration  responses  at 
70X  wing  span  due  to  excitation  by  the  trailing  edge  control  surface.  The  synthesized  approximations 
of  these  two  curves  are  very  accurate  and  practically  coincide  with  the  computed  curves.  Further  exam- 
ination of  figure  3 reveals  that  the  vertical  acceleration  response  contains  substantial  contributions 
from  modes  adjacent  to  the  flutter  mode  characterized  by  the  sharp  peak  at  13  Hz.  The  same  general  ob- 
servation was  made  for  all  three  configurations.  When  designing  a flutter  suppression  system  it  was 
found  Important  to  maximize  the  relative  response  In  the  flutter  mode  so  that  the  active  system  at  nominal 
gain  values  performs  without  undue  excitation  of  the  adjacent  modes. 


FIGURE  1.  WING/STORE  MODEL  WITH  ACTIVE  FLUTTER  CONTROL 
CONFIGURATION  (A) 


The  block  diagram  presented  in  figure  4 Illustrates  the  flutter  suppression  system  used  in  the  test 
program.  In  order  to  ensure  stability  of  the  adjacent  modes,  several  filters  are  used  to  condition  the 
feedback  signal  (V2  - Vj),  which  is  the  differential  voltage  delivered  by  two  accelerometers  located  at 
701  of  the  wing  span.  A differential  acceleration  of  one  g (9.8  m/s^)  provides  a voltage  of  approximately 
0.2V.  Two  notch  filters  at  90  Hz  and  34  Hz  were  required  to  eliminate  ground  resonances  in  the  structural 
panels  adjoining  the  actuator.  A second-order  low-pass  filter  (breaking  at  42  Hz)  was  found  to  eliminate 
resonance  in  high-order  vibration  modes  such  as  the  tip  missile  bending  mode.  The  purpose  of  the  high- 
pass  filter  is  to  minimize  the  coupling  between  the  structural  vibration  modes  and  the  rigid-body  modes 
Introduced  by  the  suspension  system.  In  the  present  program,  emphasis  was  placed  on  simplicity  In  design 
and  fabrication  of  the  flutter  suppression  system.  To  satisfy  this  requirement  It  was  found  that  the 
most  desirable  filter  for  compensation  of  the  flutter  mode  Is  a first-order  variable  gain/phase  network 
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as  shown  In  the  block  diagram.  A dial  on  the  control  panel  is  used  to  adjust  the  gain  values  In  the 
range  0. 0-2.0.  Similarly,  using  another  dial,  the  time  constant  l of  the  phase  compensat ion  can 
assume  values  between  0.0  and  0.05  seconds. 


FIGURE  4.  BLOCK  DIAGRAM  OF  THE  FLUTTER  SUPPRESSION  LOOP 

Extensive  analyses  were  performed  to  determine  the  optimal  gain-phase  settings  for  various  dynamic 
pressures  and  Mach  numbers.  The  approximations  of  the  aeroelastlc  response  at  various  speeds  and  Mach 
numbers  were  used  to  obtain  root  locus  plots  with  the  loop  gain  as  parameter.  Figure  5 presents  a 
typical  root  locus  plot  Cor  configuration  IB)  with  leading  edge  control,  at  a speed  sllghtlv  exceeding 
the  flutter  speed.  It  is  observed  that  all  poles  associated  with  structural  vibration  modes  are  located 
close  to  the  imaginary  axis.  Another  finding  is  that  the  flutter  suppression  system  has  an  adverse 
effect  on  the  first  wing  bending  mode  which  becomes  unstable  at  high  loop  gains.  Due  to  this  limiting 
factor  on  the  allowable  loop  gain,  the  flutter  suppression  system  is  able  to  stabilize  the  system  onlv 
up  to  moderate  dynamic  pressures  as  demonstrated  in  Figures  6 and  7 which  present  magnified  views  of  the 
roots  associated  with  the  first  two  structural  modes.  All  roots  associated  with  the  actuator  servo  and 
the  filters  are  sufficiently  stable  for  normal  gain  values. 

Similar  analyses  were  performed  for  the  other  configurations  with  both  leading  and  trailing  edge 
controls  at  various  Mach  numbers.  The  trailing  edge  system  for  configuration  (A)  was  found  to  be  very 
effective  in  suppressing  the  "hump"  mode  flutter.  For  configuration  (C) , as  for  (B) , the  analysis 
predicted  that  the  leading  edge  surface  would  be  more  effective  in  controlling  the  flutter  mode.  For 
both  conf lgurat ions,  however,  the  allowable  gain  was  limited  bv  the  destabilizing  effect  predicted  for 
the  root  associated  with  fundamental  wing  bending. 

4 . GROUND  VIBRATION  TEST 

The  initial  ground  vibration  test  (GVT)  was  conducted  at  Wright-Patterson  Air  Force  Base  in  Ohio 
during  April  1977.  Two  different  types  of  vibration  tests  were  performed.  The  conventional  type 
utilized  discrete,  sinusoidal  shaker  input  at  each  natural  frequency.  Relative  accelerations  were 
measured  by  a roving  accelerometer . Due  to  the  time-consuming  nature  of  this  test,  it  was  completed 
for  configuration  (A)  only.  The  second  type  of  GVT  utilized  a Hewlett-Packard  Fourier  Analyzer,  System 
5451B  operated  by  the  Dynamics  Test  Group  at  AFFDL.  The  wind  tunnel  model  was  excited  bv  random  noise 
input  to  a shaker  located  in  the  fuselage  nose.  The  response  was  measured  consecut ivel v at  preselected 
stations  and  transmitted  to  the  computer.  Almost  simultaneously  with  the  data  transmittal,  the  transfer 
function  for  each  test  point  was  computed  and  stored  on  disc.  Sampling  time  per  point  was  approximately 
120  seconds  and  the  resulting  frequency  resolution  was  0.2  Hz.  The  transfer  function  data  stored  on 
disc  were  used  to  determine  the  modal  deflections  at  approximately  85  points  for  each  of  the  three 
configurations.  To  extract  accurate  measurements  of  the  modal  damping,  a zoom  technique  was  employed 
with  data  sampled  at  a few  locations  during  approximately  bOO  seconds.  Figure  8 presents  a tvpical 
transfer  function  computed  by  the  Fourier  Analyzer.  Even  though  all  modes  generated  by  the  Fourier 
Analvzer  betwee  i 0 and  50  Hz  are  clearly  identifiable,  the  shapes  of  the  low-frequency  modes  show  poor 
correlation  with  the  corresponding  modes  measured  in  the  conventional  GVT  or  the  ones  obtained  from 
analysis.  The  results  for  higher-order  modes  are  satisfactory.  Later,  during  each  wind  tunnel  entry, 
the  conventional  vibration  test  was  repeated  to  verify  the  structural  integrity  of  the  model  and  deter- 
mine the  influence  of  the  suspension  system. 


FIGURE  7.  FIRST  BENDING  MODE  ROOT  LOCUS,  CONFIGURATION  (B)  , 
LEADING  EDGE  CONTROL  M - 0.8 
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FIGURE  8.  ACCELERATION  AT  LEADING  EDGE  PICKUP  DUE  TO  RANDOM  INPUT 
AT  FUSEUGE  NOSE,  CONFIGURATION  (A> 

Functional  checks  of  the  actuator  servos  ami  other  components  In  the  control  system  were  made  both 
In  a bench  test  and  with  the  system  installed  In  the  model.  Open-loop  tests  were  performed  to  measure 
the  servo  output  versus  input  command . The  gains  of  the  actuation  loop  were  set  for  optimum  closed-loop 
response  * which  then  was  modified  bv  various  filters  to  match  a tvplcal  aileron  actuator.  After  completion 
of  the  functional  tests,  the  flutter  suppression  loop  was  closed  to  determine  the  loop  gain  at  ground 
resonance.  Due  to  the  potential  danger  to  the  model,  precautions  were  taken  to  limit  the  hydraulic  flow 
rate  provided  bv  the  servo-valve.  Notch  filters  were  inserted  in  the  flutter  suppression  loop  to  eliminate 
critical  modes.  To  further  increase  the  available  loop  gain,  the  notch  filters  were  modified  and  a second- 
order  low-pass  filter  Inserted  during  the  second  wind  tunnel  entry. 

In  addition,  a frequency  response  test  was  performed.  The  objective  was  to  measure  acceleration 
responses  due  to  control  surface  inputs  and  store  the  data  on  magnetic  tape  for  future  processing.  From 
the  input  and  output  signals  of  the  actuator  servo,  a transfer  function  as  shown  in  figure  ^ was  derived 
and  compared  to  the  discrete  data  measured  in  the  functional  test.  The  low-f requencv  response  of  the 
actuator  servo  is  referenced  to  1.8  degrees/volt.  After  the  first  wind  tunnel  entry,  the  servo  was 
completely  redesigned  and  the  frequency  response  was  modified  to  some  extent  in  the  high-frequency  band. 

S.  WIND  TUNNEL  TESTS 

5 . 1 FI rat  Entry 

The  wind  tunnel  test  program  was  conducted  at  the  NASA/t.angley  Research  Center  in  three  separate 
entries  during  June,  August  and  December  1*177.  In  the  first  entry,  the  leading  edge  surface  was  used  for 
sweep-excitation  of  conf igurat ion  (A).  A rapidly  diverging  oscillation  of  the  leading  edge  servo  was  en- 
countered well  below  the  predicted  flutter  speed. 

In  reviewing  the  test  results  from  the  first  entry,  Northrop  decided  to  redesign  the  actuator  servos. 

A detailed  analysis^  of  the  leading  edge  servo  was  performed  considering  the  dost abi l ir tng  effect  of  the 
airloads.  Figure  10  demonstrates  that  the  root  associated  with  an  Integrator  in  the  forward  servo-loop 
couples  with  the  root  associated  with  the  airload  on  the  leading  edge  control.  With  maximum  attainable 
loop  gain,  the  analysis  shows  that  the  system  was  only  marginally  stable.  Rased  on  this  analysts,  the 
first  step  was  to  eliminate  the  integrator  and  redesign  the  servo.  At  the  same  time,  a number  of 
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FIGURE  10.  ROOT  LOCUS  OF  LEADING  EDGE  SERVO  IN  FIRST  TUNNEL  ENTRY 


precautionary  measures  was  undertaken  to  lncreaee  the  dynamic  stiffness  of  the  aervoa.  Hie  preaaurr  feed- 
back uaed  In  combination  with  poaltlon  feedback  In  the  aervo  loop  waa  eliminated  and  the  hydraulic  aup- 
ply  preaaure  waa  doubled  to  13800  kPa. 

5.2  Second  Entry 

In  the  aecond  wind  tunnel  entry,  the  flutter  boundary  of  configuration  (A)  waa  explored.  ibe  analy- 
ala  of  thla  configuration  had  predicted  a moderately  unatable  "hump"  mode  with  algnlflcant  participation 
of  the  fuaelage  first  bending  mode.  Unexpectedly  high  damping  In  the  model  support  system,  however, 
changed  the  dynamic  coupling  and  made  the  predicted  flutter  mode  stable.  Nevertheless,  the  test 
demonstrated  that  the  control  system  Increased  the  damping  and  decoupled  the  modes  part lc lpat lng  In  the 
expected  flutter  mechanism. 

During  testing  of  configuration  (B),  sweep  excitation  using  the  trailing  edge  control  surface  was 
performed  to  measure  damping  levels  and  define  the  flutter  boundary  with  the  flutter  suppression  system 
Inactive.  For  Mach  number  0.8,  a violent  wing  torsional  oscillation  was  encountered  at  a dynamic  pressure 
of  3.64  kPa.  The  model  experienced  some  damage  and  the  test  was  terminated  for  repair. 

5.3  Third  Entry 


The  third  tunnel  entry  was  Initiated  with  testing  of  a modified  configuration  (A)  having  a maaa  of 
0.2  kg  added  to  the  AIM-9E  missile  tip.  During  this  teat,  moat  damping  trends  were  obtained  using  the 
peak-hold  spectrum  method,  which  provides  Fourier-transformed  measurements  of  the  response,  filtered 
through  250  narrow-band  circuits.  The  peak  response  within  each  interval  la  reglatered  on  a screen, 
permitting  the  test  engineer  to  determine  when  the  resulting  spectrum  has  converged  and  the  data  sampling 
can  be  discontinued.  The  damping  of  a resonating  mode  Is  proportional  to  the  inverse  of  the  peak- 
hold  amplitude.  Figure  11  Illustrates  the  model  response  In  terms  of  peak-hold  spectra  of  the  wing 
root  torsion  moment  due  to  either  tunnel  turbulence  alone  or  tunnel  turbulence  In  combination  with  sweep 
excitation  by  a leading  edge  control  surface.  Although  the  flutter  suppression  system  was  not  activated 
for  the  test  points  covered  In  this  graph,  the  data  gives  an  estimate  of  the  relative  tunnel  turbulence 
and  the  threshold  amplitude  at  which  the  flutter  suppression  system  can  be  expected  to  affect  the  response. 
Figure  12  presents  a peak-hold  damping  trend  for  configuration  (A)  with  and  without  flutter  suppression. 

As  In  the  previous  tunnel  entry.  It  was  demonstrated  that  the  trailing  edge  control  system  provided  a 
significant  amount  of  damping,  but  friction  In  the  model  suspension  system  was  sufficient  to  prevent 
flutter  for  this  configuration. 

The  test  continued  with  configurations  (B)  and  tC) , which  were  designed  specifically  to  have 
severe  flutter  on-set.  Figures  13  and  14  present  measured  damping  trends  for  configuration  (B)  at  Mach 
numbers  0.6  and  0.8,  respectively.  The  model  sustained  slight  damage  at  M * 0.8  where  a "hard"  flutter 
point  was  obtained.  Subsequently,  considerable  Increases  In  the  flutter  speeds  were  achieved  using  the 
leading  edge  control  surface.  For  configuration  (B)  at  Mach  number  0.6,  a demonstrated  improvement  of 
1ST  and  a projected  Improvement  of  29X  In  the  dynamic  pressure  were  accomplished.  As  predicted  by  analy- 
sis, the  Increased  damping  In  the  flutter  mode  was  traded  for  a destabilizing  effect  on  the  first  wing 
bending  mode.  For  this  reason,  the  nominal  loop  gain  In  the  flutter  suppression  system  was  limited  to 
0.7.  The  frequency  shifts  predicted  by  the  analysis  for  system-on,  were  also  substantiated  by  the  test 
showing  a sudden  drop  In  the  torsion  frequency  as  the  dynamic  pressure  approached  Its  critical  value. 

Table  2 presents  estimated  damping  ratios  for  the  data  points  shown  In  Figure  14.  These  results  were 
obtained  by  post-test  analysis  of  analog  tape  records  using  a movlng-block/randomdec  technique  as 
described  In  reference  12. 


Dynamic 

Pressure 

(kPa) 


Damping  Ratio 


Syates>-Off 


System-tin 


Table  2.  Damping  Ratios  of  Flutter  Modes  for  Configuration  (B) 

M ■ 0.80,  Leading  Edge  Control 

Similar  teat  results  were  outalned  for  configuration  (C).  As  for  configuration  (B),  the  critical 
mode  with  the  flutter  suppression  system  on  was  the  first  wing  bending  mode.  The  effect  of  the  leading 
edge  system  on  the  peak-hold  response  Is  presented  In  Figure  15.  Since  In  thla  case  the  projected  dynamic 
pressure  at  flutter  was  5.10  kPa,  the  upper  two  plots  demonstrate  the  effectiveness  of  the  closed-loop 
system.  Similar  data  were  obtained  with  the  trailing  edge  system  In  operation.  Damping  trends  at  Mach 
0.6  are  presented  In  Figures  16  and  17  with  the  leading  and  trailing  edge  systems  on  and  off. 

Future  flutter  suppression  systems  moat  likely  would  be  tied  to  special  purpose  computers  responding 
In  an  adaptive  manner  to  counteract  the  structural  response.  An  alternative  and  technically  simpler  approach 
would  be  to  depend  on  the  flight  data  computer  to  provide  pre-programmed  control  laws  for  various  config- 
urations and  flight  conditions.  In  the  present  program,  It  was  observed  from  both  test  data  and  analytical 
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predictions  that,  as  a general  rule,  the  flutter  suppression  loop  required  increasing  phase  lag  with 
increasing  Mach  number.  The  analytls  Indicated  that,  in  some  cases,  scheduling  of  the  loop  gain  would 
be  beneficial  to  maximize  system  damping  at  the  subcritical  dynamic  pressures  and  to  stabilize  structural 
modes  adjacent  to  the  flutter  mode.  Further  testing  is  required  to  Investigate  transonic  effects  and 
improvements  of  the  control  laws. 

6.  ANALYTICAL  CORRELATION  WITH  TEST  DATA 

6. 1 Analysis  without  AFSS 

The  Air  Force  Flight  Dynamics  Laboratory  conducted  flutter  analyses  for  each  of  the  unaugmented 
configurations  using  the  FOP  module  of  the  FASTOP  (Flutter  and  Strength  Optimization  Program)  computer 
program^).  For  configuration  (A)  at  M»0.8,  flutter  was  predicted  at  a dynamic  pressure  of  6.42  kPa 
using  computed  vibration  modes  and  at  a dynamic  pressure  of  6.56  kPa  using  ground  vibration  test  data. 

It  was  noted  earlier  that  this  configuration  would  not  flutter  in  the  tunnel  due  to  unexpectedly  high 
damping  in  the  support  system.  For  Conf igurat ion  (B)  the  corresponding  numbers  were  3.64  kPa  and  4.21 
kPa.  Using  computed  modes,  the  analysis  was  about  101  conservative.  With  measured  modes  it  was  about 
4%  unconservat ive.  The  analysis  for  Configuration  (C)  at  M=0.6  using  computed  vibration  modes,  pre- 
dicted a flutter  instability  at  4.93  kPa,  2Z  below  the  projected  test  point. 

6 . 2 Analysis  with  AFSS 

To  determine  the  stability  characteristics  of  the  augmented  model,  the  AFFDL  employed  a frequency 
domain  procedure  with  a modified  Nyquist  criterion.  Figure  18  presents  a series  of  Nyquist  plots  for 
configuration  (B)  with  leading  edge  control  at  M“0.8.  At  low  dynamic  pressure  where  the  response  is 
small,  the  loop  corresponding  to  the  critical  mode  is  small  and  located  to  the  right  of  the  origin.  As 
the  flutter  speed  is  approached,  the  Nyquist  loop  gets  larger  and  approaches  a curve  of  infinite  radius 
that  turns  in  the  clockwise  direction.  At  flutter,  the  loop  degenerates  to  a straight  line.  As  the 
flutter  speed  is  surpassed,  the  critical  mode  encircles  the  origin  in  a counter-clockwise  manner 
indicating  that  the  model  is  stable.  For  higher  speeds  (above  flutter),  the  Nyquist  loop  becomes  smaller 
until  it  no  longer  encircles  the  origin.  At  this  point,  suppression  of  the  flutter  mode  is  lost.  For 
configuration  (B),  loss  of  flutter  control  was  predicted  at  a dynamic  pressure  of  5.27  kPa.  This 
represents  an  increase  of  31%  in  flutter  dynamic  pressure  over  the  unaugmented  case.  Experimental  re- 
sults at  this  Mach  number  indicated  a projected  increase  of  15%  using  the  peak-hold  trend  and  a projected 

increase  of  23%  using  the  est invited  damping  ratios  presented  in  Table  2. 

The  analyses  for  the  augmented  configuration  (A)  at  M*0.8  showed  that  the  critical  dynamic 

pressure  could  be  increased  by  about  33%  using  either  the  leading  edge  or  the  trailing  edge  surface. 

Since  this  configuration  was  stable  in  the  tunnel,  it  is  difficult  to  make  a meaningful  comparison. 

For  the  augmented  configuration  (C)  with  leading  edge  control  the  analysis  at  M*0.6  showed  good 
stability  characteristics  up  to  dynamic  pressures  exceeding  9.58  kPa  which  is  48%  above  the  unaugmented 
case.  The  test  results,  however,  show  a projected  improvement  of  only  17%. 

In  general,  control  surface  aerodynamic  force  and  moment  coefficients  predicted  by  theory  are  high 
when  compared  to  experimental  data.  Since  experimental  information  on  this  model  regarding  control 
surface  aerodynamics  was  not  available,  there  was  no  attempt  at  this  time  to  reduce  the  magnitude  of 
the  control  surface  aerodynamics.  Unmodified  control  surface  aerodynamics  tend  to  predict  a more 
efficient  active  flutter  suppression  system.  This  may  be  a partial  explanation  of  the  differences  ob- 
tained between  the  analysis  and  test  results. 

7.  CONCLUDING  REMARKS 

Although  many  problems  remain  to  be  solved  in  order  to  make  a flutter  control  system  truly  adapt- 
ive, the  present  test  program  has  demonstrated  that  active  suppression  of  wing/store  flutter  is  feasible 
for  practical  application.  As  far  as  known,  this  was  the  first  time  that  a leading  edge  surface  was 
used  as  the  single  active  device  for  flutter  control.  Both  the  leading  and  the  trailing  edge  sur- 
faces were  used  independently  to  suppress  a single  flutter  mode.  There  are  several  items  that  remain 
as  future  tasks.  At  present,  the  Air  Force  is  reviewing  plans  to  use  the  existing  model  for  demonstration 
of  improved  control  laws  and  various  adaptive  control  schemes  for  flutter  suppression.  Another  item 
considered  for  future  testing,  is  the  combined  application  of  leading  and  trailing  edge  controls. 
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14.  Abstract 


Air  Forces  in  many  countries  have  to  face  problems  of  aeroelasticity  and  flutter  with 
aircraft  carrying  more  and  more  stores.  The  two  papers  of  this  Report  presented  to  the  Sub- 
committee on  Aeroelasticity  of  the  Structures  and  Materials  Panel  during  the  46th  Meeting 
of  the  Panel  deal  with  two  different  aspects  of  the  problem. 

The  paper  of  A.Lotze  is  concerned  with  the  many  configurations,  symmetric  or 
asymmetric,  that  may  occur  and  with  their  consequences  on  the  natural  modes  of  the 
structure.  It  clarifies  the  difficulties  one  has  to  face  to  clear  the  (light  domain  for  all  flight 
configurations  and  proposes  useful  approaches. 

The  paper  of  C.llwang,  B.A.Winther,  T.E.Noll  and  M.G. Farmer  deals  with  the  difficult 
problem  of  flutter  suppression;  it  exposes  the  approach,  the  design,  and  the  wind  tunnel  tests 
of  the  model  of  a fighter,  carrying  stores  and  equipped  with  a flutter  suppression  device. 


The  two  papers  are  of  great  interest  to  aeroclasticians  and  may  give  useful  help  to  the 
designer. 


- The  paper  of  A.Lotze  is  concerned  with  the  many  configurations,  symmetric  or  The  paper  of  A.Lotze  is  concerned  with  the  many  configurations,  symmetric  or 
asymmetric,  that  may  occur  and  with  their  consequences  on  the  natural  modes  of  the  asymmetric,  that  may  occur  and  with  their  consequences  on  the  natural  modes  of  the 
structure.  It  clarifies  the  difficulties  one  has  to  face  to  clear  the  flight  domain  for  all  structure.  It  clarifies  the  difficulties  one  has  to  face  to  clear  the  flight  domain  for  all 
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